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Abstract: We construct a new, simple phenomenological model along the lines of 
AdS/QCD. The essential new ingredient is the brane-antibrane effective action includ- 
ing the open string tachyon proposed by Sen. Chiral symmetry breaking happens because 
of tachyon dynamics. We fit a large number of low-spin meson masses at the 10%-15% 
level. The only free parameters involved in the fits correspond to the overall QCD-scale 
and the quark masses. Several aspects of previous models are qualitatively improved. 
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1. Introduction 

Understanding the strong dynamics underlying many observations related to the strong 
interaction remains an unsolved problem. Much progress has been made to date with 
different methods but new insights are always welcome. A recent development that has 
led to reconsidering the strong interaction has been the AdS/CFT duality. This has been 
applied to obtain new insights on QCD phenomenology. In this paper, we focus on the 
meson spectrum - (see [|| for a review of the gauge-gravity literature on the issue). 

There are two main ways to address the problem at hand. Top-down approaches use 
string theory from first principles in order to build dual theories as close as possible to 
QCD. Notable examples are 0], 0. On the other hand, bottom-up approaches, starting 
from the works |5|, ^, use known QCD features to develop holographic models that are 
only inspired by string theory. The model we present here is more of the second kind, but 
includes the main stringy ingredients we expect from first principles, namely the effective 
action controlling the chiral dynamics. Our main observation is that merging this stringy 
input of top-down holographic models for flavor allows us to improve the existing bottom- 
up models both qualitatively and quantitatively. 

2. The model 
Background and action 

In 1^, it was shown, quite generally, that effective actions for brane-antibrane systems 
derived from string theory |^] encode a set of qualitative features related to chiral symmetry 
breaking and QCD at strong coupling. The goal of this paper is to build a concrete model 
within the framework of Q . We will consider the simplest smooth gravitational background 
that is asymptotically AdS, while having a confining IR in the same spirit as Q. This 
turns out to be the AdS^ soliton, which was shown to be a solution to the two-derivative 
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approximation of subcritical string theory and used as a toy model for certain aspects of 
4d Yang-Mills in |10]. The metric reads: 



dsl = ^ [dxl^ + f^^^dz^ + /a V] (2.1) 

with /a = 1 — The coordinate rj is periodically identified and z G [0, za]- The dilaton is 
constant and we do not write the RR-forms since they do not play any role in the following. 
We now consider a D4-D4 pair, located at fixed rj in this background^. We write the action 
proposed by Sen [^j as 

S = - j d^xdzV{\T\) (^V-detAL + V-detAii) (2.2) 

The objects inside the square roots are defined as: 



^MN = 9MN + W [2F];^^ + {{DmTT (DnT) + (M ^ iV))J (2.3) 

where M,N = 1, ... ,5, the field strengths = OmA^^ — d^A^^j and the covariant 

derivative of the tachyon is DmT = {du + ~ *^m)^- The active fields in (^^), ( ^^ ) 
are two 5-d gauge fields and a complex scalar T = re*^, which are dual to the low- lying 
quark bilinear operators that correspond to states with J^'-' = 1 ,1"'""'",0 '",0^"'' -; see 
for details. In the action of Q, the transverse scalars (namely rj in the present case) 
are also present. We have discarded them when writing ( |2.3| ) since they do not have any 
interpretation in terms of QCD fields. Accordingly, even if the background ( p.l| ) is six- 
dimensional, the holographic model for the hadrons is effectively five-dimensional and, in 
fact, its field content coincides with those of IH, |6|. For the tachyon potential we take, as 
the simplest possibility, the one computed in boundary string field theory for an unstable 
Dp-brane in fiat space [||], although one should keep in mind that this expression for V is 
not top-down derived for the present situation. In the present conventions, V = lCe~^'^ , 
where fC is an overall constant that will play no role in the following since it does not 
enter the meson spectrum computation (it is important though, in the normalization of 
correlators when computing for instance decay constants |^], Q, |Q). The tachyon mass 
is = —-^ and we will impose: i?^ = Go;' in order to have m^R"^ = —3. This should 
not be interpreted as a modification of the background due to the branes, but just as a 
(bottom-up) choice of the string scale that controls the excitations of those branes, such 
that the bifundamental scalar T is dual to an operator of dimension 3, as in [Q. Since 
the AdS radius is not parametrically larger than a' , the two-derivative action cannot be a 
controlled low enegy approximation to string theory. This is the main reason why a model 
of this kind cannot be considered of top-down nature. Notice the value of we take differs 
from the one used in [fi!o|l. 



.(0 



In five-dimensional holographic models of QCD, the flavor branes are expected to be a D4-D4 system, 
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The tachyon vacuum and chiral symmetry breaking 

As shown in 0, an essential ingredient of the present framework is that the generation of 
the correct flavor anomaly on the flavor branes requires the tachyon modulus r to diverge 
somewhere. Therefore, r must have a nontrivial vev which breaks the chiral symmetry. 
From the action ( ^.2[) we obtain the equation determining t{z): 



where the prime stands for derivative with respect to z. Near z = 0, the solution can be 
expanded in terms of two integration constants as: 

T = ClZ + ^cfz^ log Z + C3Z^ + 0{z^) (2.5) 

b 

where, on general AdS/CFT grounds, ci and C3 are related to the quark mass and con- 
densate (see for a careful treatment). From ( ^.4] ), we find that r can diverge only at 
z = z/y. There is a one-parameter family of diverging solutions in the IR: 



C 13 



3 



3 67rC^^^"^^'^ + --- ^^-^^ 



{ZA - z) 20 

The interpretation is the following: for a given ci (namely quark mass^) fixed in the UV 
(near z = 0), the value of C3 (namely the quark condensate) is determined dynamically 
by requiring that the numerical integration of (|2.4D leads to the physical IR (near z = za) 
behavior (p.6|). Hence, for any value^ of ci, one can obtain numerically the function for the 
vev (r). 

Meson spectrum: Numerical results 

There is a rather standard method for computing the meson spectrum in holographic 
models, see M for a review. Each bulk field is dual to a boundary operator and its 



linearized perturbation can be obtained after expanding (2^). By looking at normalizable 
fluctuations of the bulk fields, one typically encounters discrete towers of masses for the 
physical states with the corresponding quantum numbers. Thus, for a fixed value of ci, 
we have a Sturm-Liouville problem for each bulk mode. This can be solved numerically, 
using a standard shooting technique. By computing the different towers at different values 
of ci, we found the following expressions to be very good approximations to the numerical 
results, in the range < ci < 1 where we could perform the numerics reliably. For the 
vectors: 

ZA m^y^ = 1.45 + 0.718C1 , za = 2.64 + 0.594ci 
ZA mP = 3.45 + 0.581C1 , za m^y'^ = 4.13 + 0.578ci 
ZAmP = 4:.72 + 0.577ci , za mj^^ = 5.25 + 0.576ci. (2.7) 



^For the present work, we will just use that ci is proportional to m,. Finding the proportionality 
coefficient requires normalizing the action and fields as in |Q. 

^In practice, we have been able to perform numerics in a controlled manner only for < ci < 1. 
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For the axial vectors: 



= 1.93 + 1.23ci , za mj^ = 3.28 + 1.04ci 
ZAmJ^ = 4.29 + 0.997ci , za^J^ = 5.13 + 0.975ci 
Zf, m^^ = 5.88 + 0.962C1 , zk m^^ = 6.55 + 0.954ci . (2.8) 



For the pseudoscalars: 



zi^m^p^ = \jlAlc{ + 5.32ci , ZKm>p' = 2.79 + 1.16ci 
zi^ mP = 3.87 + 1.08ci , za = 4.77 + 1.04ci 

= 5.54 + l.Olci , ZAmi^^ = 6.24 + 0.997C1. (2.9) 



For the scalars: 



ZA m^p = 2.47 + 0.683C1 , za m^p = 3.73 + 0.488ci 
z^ = 4.41 + 0.507ci , za mlf^ = 4.99 + 0.519ci 
z^mjf^ = 5.50 + 0.536C1 , za mjf^ = 5.98 + 0.543ci. (2.10) 

It turns out that meson masses increase hnearly with ci . Namely, they increase linearly with 
the bare quark mass, as expected from an expansion in niq and in qualitative agreement 
with lattice results, see for instance [15|,[^],[|l^]. The exception, of course, is the first 



pseudoscalar for which itItj is proportional to ^/m^ (for small ruq), as expected from the 
Gell-Mann-Oakes-Renner relation. Actually, the behaviour 171^^ = ^/bniq + drn^ was also 



found in the lattice [15|. 



3. Fitting the meson spectrum 

We now proceed to make a phenomenological comparison of the results of (^.7])- (|2.10 ) to 



the experimental values quoted by the Particle Data Group [18|. Obviously, we can only 
model those mesons with J^"^ = l~,l++,0-+,0++. From |l|], we will just take the 
central value quoted for each resonance. We do not discuss decay widths here (in the strict 
Nc ^ oo limit they are of course zero). 

Isospin 1 mesons 

We start by looking at mesons composed of the light quarks u and d. In particular, we 
discuss the isovectors. In table |^, we show all the mesons listed in the meson summary 



table of 1 18 1 under light unflavored mesons which have isospin 1 and the J 's present in 



our model. The only exception is ao(980), which is considered to be a four-quark state |18| 



We have fitted the parameters of the model to these observables by minimizing the rms 



1 



error e^ms = 100 x {'§') j ' where n = 8 — 2 = 6is the number of observables 

minus the number of parameters. We obtain for the parameters 



-1 



522MeV , ci,i = 0.0125 (3.1) 
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Meson 


IVleasureu (IVleVj 


Modei (MeVj 


1— 


p(770) 


775 


762 


p(145Uj 


14d5 


1379 




i7z(J 


loOb 


1++ 


ai(1260) 


1230 


1015 


0-+ 




135.0 


135.1 


7r(1300) 


1300 


1462 


7r(1800) 


1816 


2026 


0++ 


ao(1450) 


1474 


1295 



Table 1: A comparison of the results of the model to the experimental values for light unflavored 
meson masses. 



with Erms = 12%. 

In table we display the resonance masses with isospin 1 listed in [ p!^ ] under other 
light unflavored mesons. These are namely states considered as "poorly established that 
thus require confirmation" . For the results given by our model, we use ( |3.1D and therefore 
no further parameter is fitted here. For this set of observables we get e^ms = 24%, where 
we have inserted n = ll — = 8. One should keep in mind that it is plausible that some 
of these "unconfirmed" states may not be real or may be misinterpreted as part of the 
meson towers. In this sense, our model seems to favor the /o(2150) as the fourth member 
of the /3-meson tower ^. We have not included p(1570) in table § because its excitation 
number is smaller than /9(1700), which was included in the previous fit. In case /9(1570) 
gets confirmed as a member of this tower, the fit should be redone. We observe that the 
model tends to consistently overestimate the masses of the excited axial vectors and pions. 
This is connected to the fact the model yields a Regge slope for axial mesons larger than 
the one for the vectors If the experimental results of table ^ are confirmed, one should 
think of improving the model in order to ameliorate this discrepancy. 



ss states 

A nice feature of the present model is that it incorporates the dependence of the hadron 
masses on the quark mass. This allows us to study ss states. More precisely, it allows us to 
discuss "hypothetical states" with quark content ss assuming no mixing with other states. 
In the real world, the mixing for pseudoscalars and scalars is important (see chapter 14 of 
flSll), and therefore it is not possible to compare directly the outcome of the model to the 



experimental results. Nevertheless, as in [19], we can estimate the masses of these "hypo- 



thetical" ss mesons from the light-strange and light-light mesons. Then, using quotation 
marks for the hypothetical states, and using the quark model classification (table 14.2 of 



11]) we can write m("?7") = yj2m]^-ml, m("(/)(1020)") = 2m(K*(892)) - m(p(770)), 
m("r/(1475)") = 2m(E:(1460)) - m(7r(1300)), etc. Keeping the value of found in (U), 



''indeed, there is much more experimental evidence for p(2150) than for p(1900) or p(1570). We thank 
S. Eydelman for very useful explanations. 
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Meson 


iVleasurea (MeVj 


Modei (MeVj 


1— 




19UU 


2159 


p(215Uj 


215U 


2467 


p(2270j 


2270 


2746 


1++ 


ai(lo4Uj 


1647 


1721 






2245 


«! (zuyb j 


zUyb 


zbob 


ai(2270) 


2270 


3073 


ai(2340) 


2340 


3423 




7r(2070) 


2070 


2493 


7r(2360) 


2360 


2899 


0++ 


ao(2020) 


2025 


1952 



Table 2: A comparison of the results of the model to the experimental values for other light 
unflavored meson masses. 

we fit the value of ci associated with the strange quark to the "experimental" values of 
table ^, obtaining. 

ci,, = 0.317 (3.2) 
The rms error for this set of observables (n = 6 — 1) is £rms = 10%. 





Meson 


Measured (MeV) 


Model (MeV) 


1— 


"(/)(1020)" 


1009 


876 


"(/)(1680)" 


1363 


1474 


1++ 


7i(1420)" 


1440 


1210 


0-+ 


"?7" 


691 


725 


"ry(1475)" 


1620 


1647 


0++ 


"/o(1710)" 


1386 


1403 



Table 3: A comparison of the results of the model to the hypothetical states with ss quark content. 



FYirther comments 

We end this section by commenting on the estimates given by the model on two other 
physical quantities. Without fixing the proportionality coefficient between ci and the 
quark mass, from ( p. 11 ), (p.2|), we can infer the ratio of the strange quark to the light quark 
mass: ^^"^^ ~ ~ 25. Moreover, the background studied here experiences a first order 
deconfinement phase transition in full analogy with The deconfinement temperature 
is given by Tdeconf = ^^^k ^ MeV. Both of these values are close to the experimental 
values. 
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4. Conclusions 



We have built a new phenomenological model for the meson sector of QCD. In this paper 
we have discussed the mass spectrum. We note the simplicity of the construction, whose 
essential point is the use of Sen's action [0] including the open string tachyon field. We 



have applied it to one of the simplest backgrounds exhibiting confinement |10|. Despite 



the minimal input, we have found the following interesting qualitative properties: 

• The model includes towers of excitations with J^^ = 1^^, 1^^, ^, 0^^, namely all 
low-lying operators that do not need a dual excited stringy state. 

• Chiral symmetry breaking is consistently realized. Moreover, the value of the quark 
condensate is computed dynamically and is not a tunable input. Hence, the num- 
ber of tunable parameters coincides with those present in QCD: they are just the 
dynamically generated scale and the quark masses. 



We find Regge trajectories for the excited states oc n, as in the soft wall model 
[pO|]. This allows good predictions for the higher excitations, as opposed to the hard 
wall model |5|, ^. Notwithstanding, the Regge slope for axial vectors is bigger than 
the one for vectors. This fact requires further study. 

• Our model incorporates the increase of the vector meson masses due to the increase 
of quark masses, as rrip ~ fci + /c2 for small rriT^. 

Previous AdS/QCD models present some of these properties, but as far as we know, 



no existing model is able to encompass all of them, see |21] for recent related discussions. 
We briefly comment on the three benchmark models: the Sakai-Sugimoto model Q misses 
the first and third points listed above, the hard wall model [S, 0] misses the third one and 



partially the second one; and the soft wall model [20| misses the second one. All of these 
models [^, ^, ^, ^ and variations thereof fail to get the fourth point (although it is worth 
mentioning that D3D7 models with abelian flavor symmetry do capture the physics of this 
fourth point, see section 6.2.3 of 

Moreover, the quantitative matching shown in tables || and ^ with the central values 
of the meson resonances is excellent, at the 10%-15% level. This is a typical accuracy of 
AdS-QCD-like models (a recent example is [^], which accounts for excited spin states of 
the p and uj families). Since the systematic error produced by quenching is of the order 
of 10% |2^] and the differences between quenched lattice computations with iVc = 3 and 
Nc = oo are again of the order of 10% |16, it would be unexpected to get a better 



accuracy from any model of the kind presented here. 

It would be of utmost interest to generalize the set-up to the non-abelian case, allowing 
several quark flavors, but this is beyond the scope of the present work. 
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